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Abstract. The nonequilibrium stationary state of an irreversible spherical model 
is investigated on hypercubic lattices. The model is defined by Langevin equations 
similar to the reversible case, but with asymmetric transition rates. In spite of being 
irreversible, we have succeeded in finding an explicit form for the stationary probability 
distribution, which turns out to be of the Boltzmann-Gibbs type. This enables one to 
evaluate the exact form of the entropy production rate at the stationary state, which 
is non-zero if the dynamical rules of the transition rates are asymmetric. 



PACS numbers: 05.10.Gg, 05.70.Ln, 75.10.Hk 



Irreversible spherical model and entropy production rate 



2 



1. Introduction 

The spherical model was introduced by Kac pLj as a modification of the Ising model 
in which discrete spin variables are replaced by continuous ones, but subjected to the 
spherical constraint, which is a condition that ensures the thermodynamic properties 
of this system for any temperature. The critical behaviour of this model was first 
analysed by Berlin and Kac [2], and the exact solution can be found in any dimension d. 
The model displays a continuous phase transition for d > 2, with non-classical critical 
behaviour for 2 < c? < 4 and mean field properties for d > 4. The rich critical behaviour 
[3], together with the establishment of many exact results, has made the spherical model 
a nice laboratory for statistical mechanics methods. 

Being defined in a static way by the Boltzmann-Gibbs probability distribution, the 
spherical model has no dynamics. However, a dynamics can be assigned to the model 
by the introduction of a set of Langevin equations [1], which will rule the time evolution 
of the spin variables now transformed into stochastic variables. The Langevin equations 
have additive white noise and the deterministic parts are linear in the stochastic variables 
[5], [6] . These variables, as in the static case, are associated to the sites of a regular lattice 
and, in addition, they are subject to the spherical constraint. The stationary probability 
distribution of the associated Fokker-Planck equation turns out to be the Boltzmann- 
Gibbs probability distribution of spherical model. 

The time dependent behaviour of the dynamic spherical model defined by the 
Langevin equation has been examined in a series of papers [3, El EJ. In these works, 
the relaxation to the thermodynamic equilibrium was investigated through two-point 
functions (autocorrelation and response function), which enables one to quantify a 
distance of the system from the equilibirum state [10]. This approach is based on an 
extension of the fiuctuation-dissipation theorem to non-equilibrium states [HI [12] . We 
remark that the nonequilibrium situations analysed in these papers, in which the system 
is relaxing to the equilibirum, should be distinguished from the ones that concern us 
here, namely, the situation in which the system finds itself in a nonequilibrium stationary 
state. 

The deterministic part of the Langevin equations, which we call force, may be 
understood as the gradient of the Hamiltonian defining the spherical model. In other 
words, the force is conservative. Since the Hamiltonian is a quadratic form, the force 
is linear in the stochastic variables so that the linear coefficients make up a symmetric 
matrix. In this paper, we consider Langevin equations for which the forces are still linear, 
but the coefficients lose the symmetric property becoming no conservative. The set of 
Langevin equations with these nonconservative linear forces, together with the spherical 
constraint, defines the irreversible spherical model. We show here, that in spite of the 
irreversibility, the stationary probability distribution can be written as being of the 
Boltzmann-Gibbs type. The description of the stationary distribution by a Boltzmann- 
Gibbs type function has already been found in models with Ising spin variables that 
lack detailed balance [131 HI [161 [H] . 
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In the stationary state, the system is no longer in the state of thermodynamic 
equihbrium because the forces are nonconservative. In this case, there will be a 
continuous production of entropy. The second purpose of this paper is to calculate the 
production of entropy in the stationary state, which, as we shall see, can be done exactly. 
The entropy production rate for systems described by a set of Langevin equations, or by 
the associated Fokker-Planck equation, can be obtained from an expression introduced 
by Tome [18], and also considered by van den Broeck [19], which was derived from an 
expression advanced by Schnakenberg [20] for systems described by a master equation. 
The critical behavior of the entropy production rate is shown to be similar to that of 
the energy of the equilibrium spherical model. 



2. Spherical model 

The spherical model [21 [3] is defined as follows. On a d- dimensional hypercubic lattice, 
with sites and periodic boundary conditions, a continuous spin variable cxr is attached 
to each site r of the lattice. The usual nearest neighbour interaction Hamiltonian is 
written as 

'H{a) = -^^Je(rr(7v+e + fJ-^fTr, (l) 
re r 

where the summation in e is over the d orthogonal unit vectors that define the d- 
dimensional hypercubic lattice. In a cubic lattice, for instance, these unit vectors are 
ei = (1,0,0), 62 = (0,1,0) and = (0,0,1). We are considering the anisotropic case 
in which the interactions are distinct for distinct directions. The symbol a stands for 
the set of configurations {ar} of the spins and Je and /i are parameters. 
The probability distribution of configuration a is given by 

Pia) = ie-^^^'^), (2) 

where /9 = 1/kBT, ks being the Boltzmann constant and T the temperature. The 
parameter /i is not free, but shoud be chosen such that 

T 

which is called (mean) spherical constraint [211 122] . 

The dynamics of the (mean) spherical model may be formulated through the 
Langevin equation 

^ = /r(^)+^r(t), (4) 

where the force /r(o") is given by 



or 



= -'e(0-r+e + (Tr-e) " SyUdr. (6) 
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As usual, the noise term rij.{t) has the properties 

M)) = and {7]rit)r]r'it')) = 2T5r,r'6it - 1'), (7) 

where F = fc^T and T is identified with the heat-bath temperature. 

The time evolution of the probability P{a,t) of state a at time t is given the 
Fokker-Planck equation 

^-E^^^(M). (s) 

r 

where Jr{a^t) is the probability current, given by 

Jr(a, t) = Ua)P{a, t) - r^P{a, t). (9) 

The probability distribution given by equation ([2]) is the stationary solution of the 
Fokker-Planck equation. In fact, in the present case, each probability current at the 
stationary state, 

Jr(a) = /r(a)P(a)-F^P(a), (10) 
vanishes, and we may say that the system is in thermodynamic equilibrium. 

3. Irreversible spherical model 

In order to induce an irreversibility, and inspired by the Langevin equation (jlj), we 
introduce the irreversible dynamics by 

^ = /r(0-)+r^r(t), (11) 

where now the forces are given by 

/r(o-) = ^(JeCTr+e + J-eO"r-e) - S/iCJr, (12) 
e 

and cannot be anymore derived from a Hamiltonian unless Je = J-e for all e. The 
parameters Jg and J_e, in this context, should be understood as the strengths of the 
transition rates of the Markovian process defined by the Langevin equation, and not as 
an exchange integral entering the Hamiltonian as in the reversible case. Notice that, as 
before, is not free but is a time dependent paramenter that should be chosen so that 
the constraint ([3]) is fulfilled. 

The Fokker-Planck equation has the same form as before, 

dP{a,t) \^ d ^. , ^ d 



/r(a)P(a,t)-F^P(a,t) 



(13) 



but now the forces /r(o") are nonconservative and given by ( |T2i) . In the stationary state, 
the probability current 

Jr{a) = Ua)P{a) - T^P{a) (14) 
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does not vanish anymore, although the stationarity condition, 

E^J.W = 0. (15) 

r 

is fulfilled for the stationary probability distribution P{cr). 

The stationary probability distribution ^(0") is obtained by assuming a form similar 
to ([2]), namely, 

P(a) = Ce*(") with ^{a) = J2Yl ^e^rf^r+e " ^ J] ^r> (16) 

re r 

where the summation is over the nearest neighbor pairs and A and {-Be} are parameters 
to be found. We start by writing the stationary Fokker-Planck equation f lT^ in the form 

E^rM=0, (17) 
r 

where 

Of d'^ ^fd^^' 

9r{(^) = + fr^ — r— - - r 



dffr dor da^ \dar 
which was obtained after dividing the stationary equation (fT5|) by -P(cr). The 
substitution of ^{cr), given by (ITB]) . and /r(o"), given by (IT^ . into (IT5]1 shows that 
g{a) is a quadratic form in the variable Ur plus a constant. This constant is 2(/x — TA), 
and should vanish. We conclude, therefore, that 

A = ^. (19) 

Using this result, gr{(y) becomes the quadratic form 

fl'r(o-) = 5e(Ce O^r+e + C-e f^r-e) " 2/1 ^ (Tr(Ce 0"r+e + C_e C^r-e) 

e e 

+ ^ 5e(Ce + C_e)0-r-e0-r+e + ^ 5e(Ce (Tr+e + C_e (Tr-e) (c^r+e' + 0-r-e')> (20) 

e e.e' 

eXe' 

where 

Ce = Je — T-Be and C_e = -'^-e " Ti^e- (21) 

The trivial solution of (fT7|) is obtained by setting Cg = and c_e = 0, which gives 
Je = TBe = J-ei leading us back to the reversible model. To get a nontrivial solution, 
we substitute the expression fl20|) into f|T7|) and rewrite it in the form 

^ ^ fie(Ce + C_e)crr - 2A ^ XI '^^'^ + C_^)ar<yv+e 
re re 

+ X X 5e(Ce + C_e)o-r-eO"r+e + X] ^'^'^^'^ + C_e)(o"r+e + 0"r-e)o"r+e' = 0, (22) 

re r e,e' 

eXe' 

which is solved by setting Cg + c_e = 0, leading to the condition 

5e = ^(Je + J-e). (23) 
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Therefore, the irreversible model defined by the equations fllip and fll2p . which 
embodies the parameters Je, J-e and /i, has a stationary state of the Boltzmann- 
Gibbs type given by ( fT6l) with the parameters {-Be} and A given by (flQl) and (l23l) . 
It is worthwhile to notice that a totally asymmetric dynamics is obtained by setting 
J_e = 0, in which case = Jef^T, a result valid in any dimension. The totally 
asymmetric dynamics has been shown to exist in systems with Ising spin variables in 
one and two dimensions [131 ■ 



4. Entropy production rate 



The variation of entropy 5 of a system with time can be splitted into two parts as 

f^n-*, (24) 

where 11 is the entropy production rate and $ is the entropy fiux from the system to 
the environment. Following [TH], the entropy production rate 11 of a non-equilibrium 
system governed by a Fokker-Planck equation can be evaluated by 

^' r^J P{a,t) ^ ' 

From the definition of entropy, 

S{t) = -A;B5^P(a,t)lnP(a,t), (26) 

cr 

we get from (HM the following expression for the entropy flux [18] : 

'^W = f E / daJM,t) frier). (27) 

r 

In the stationary state one has 11 = $, and we can use either expression (!25|) or (!27j) to 
calculate the entropy production rate. 

From now on, we restrict ourselves to the simple case in which and J_e are 
independent of e, that is, 

Je = J and J_e = J', (28) 

but J ^ J', and the parameter = B being independent of e, which leads to 
7+7' 

Using the results of the previous section in the expression ( l25l) . the entropy 
production rate per site U* = Il/N can be evaluated as 

n* = ^(go-g2)(J-Jr, (30) 

where go and q2 are defined by 

go = (cTr) and g2 = (ar-eCTr+e), (31) 
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Figure 1. Plot of 11*/ J* as a function of 6* in d = 1, = 2 and d = 3. 



with e being any one of the unit vectors. Notice that 11* vanishes when the reversibihty 
condition J = J' is satisfied, as expected. 

Using the stationary probabihty distribution ( !T6|) . we get the following results for 
go and 



10 = ^ I d^e-^Hmr 



and 



?2 = ^/ d^e-^^m)t'm, 



(32) 



(33) 



where the modified Bessel function of first kind of order n is denoted by /„. The 
parameter fj, should ensure the spherical constraint, which is simply (0"^) = 1 or go = 1- 
Setting the right-hand side of equation (l32ll equal to 1, one has B as an implicit function 
of A. If we define 9 = 1 / B = 2T / {J + J') , the production of entropy per site becomes 



where 



n* 



J* 



^(1-^2) J* 



(J - J 



/\2 



(34) 



(35) 



J+J' 

Note that the quantity 11*/ J* is a function of 6 only, and a graph 11*/ J* x ^ is 
plotted in Figure [T]for dimensions d = 1,2, 3. In one dimension, an analytical expression 
is available, and is given by 



J* 8 V / 



(36) 
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When ^ — > 0, the quantity 11* /J* aproaches a constant (in d dimensions), which is 
given by 

1 

2 



Cd 



d 



2 

1 - - 

TT 



d 



1 - COS 2ki d'^k 



d 



d>3 



(37) 



^JBzd- Eti COS fc, (27r)'^ 
where the integration above is over the Brillouin zone (BZ). For d = 3, we get 
C3 = 0.314762 Notice that in d > 3, the quantity n*/J* equals the constant q 
in the ferromagnetic phase, 9 < 9c, where 9c is given by 

d'^k 



11/" 1 

Tc'ljBzd-Y.U 



COS 



(27r) 



(38) 



For d = 3, we get = 7.913552 



5. Conclusion 



In this paper, we have investigated a system in a nonequihbrium stationary state. 
The (mean) spherical model is a suitable laboratory where many exact results are 
available, and we have succeeded in finding an exact form for the probability distribution, 
despite the fact that the system is not in equilibrium state (as testified by the non-zero 
value of the entropy production). It is worthwhile to mention that the probability 
distribution found is of the Boltzmann-Gibbs type. The knowledge of this particular 
form allowed us to explicitly evaluate the stationary entropy production. The origin of 
the nonequihbrium behaviour in our work, which is responsible for the non-zero entropy 
production, goes back to the unbalanced transition rate to the opposite direction, 
Je 7^ J-e (for any e). If, on the other hand, the condition Jg = J_e is satisfied for 
any e, the stationary entropy production vanishes. 
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